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Abstract 

<: 

■ Wideband spectrum sensing detects the unused spectrum holes for dynamic spectrum access (DSA). Too high 
sampling rate is the main problem. Compressive sensing (CS) can reconstruct sparse signal with much fewer 
randomized samples than Nyquist sampling with high probability. Since survey shows that the monitored signal 
is sparse in frequency domain, CS can deal with the sampling burden. Random samples can be obtained by the 

^ ■ analog-to-information converter. Signal recovery can be formulated as an LO norm minimization and a linear 

in" 

■ measurement fitting constraint. In DSA, the static spectrum allocation of primary radios means the bounds between 

o 

■ different types of primary radios are known in advance. To incorporate this a priori information, we divide the 

■ whole spectrum into subsections according to the spectrum allocation policy. In the new optimization model, the 
minimization of the L2 norm of each subsection is used to encourage the cluster distribution locally, while the LO 



■ norm of the L2 norms is minimized to give sparse distribution globally. Because the L0/L2 optimization is not 

convex, an iteratively re-weighted L1/L2 optimization is proposed to approximate it. Simulations demonstrate the 
proposed method outperforms others in accuracy, denoising ability, etc. 
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1. Introduction 

Cognitive radio (CR) is a very promising technology for wireless communication. Radio spectrum is 
a precious natural resource. The fixed spectrum allocation is the major way for the spectrum allocation 
now. In order to avoid interference, different wireless services are allocated with different licensed bands. 
Currently most of the available spectrum has been allocated. But the increasing wireless services, especially 
the wideband ones, call for much more spectrum access opportunities. The allocated spectrum becomes 
very crowded and spectrum scarcity comes. To deal with the spectrum scarcity problem, there are several 
ways, such as multiple-input and multiple-output (MIMO) communication [IJ, ultra-wideband (UWB) 
communication beamforming Q 01, relay Q, and so on. Investigation demonstrates that most of the 
allocated bands are in very low utility ratios CR is proposed to exploit the under-utilization of the radio 
frequency (RF) spectrum. It is a paradigm in which the cognitive transmitter changes its parameters to 
avoid interference with the licensed users. This alteration of parameters is based on the timely monitoring 
of the factors in the radio environment. 

Spectrum sensing is one of the main functions of CR. It detects the unused frequency bands, and 
then CR users can be allowed to utilize the unused primary frequency bands. Current spectrum sensing 
is performed in two steps |I3: the first step called coarse spectrum sensing is to efficiently detect the 
power spectrum density (PSD) level of primary bands; the second step, called feature detection or multi- 
dimensional sensing [I8]|, is to estimate other signal space accessible for CR, such as direction of arrival 
(DOA) estimation, spread spectrum code identification, waveform identification, etc. 

Coarse spectrum sensing requires fast and accurate power spectrum detection over a wideband and 
even ultra-wideband (UWB). One approach utilizes a bank of tunable narrowband bandpass filters. But it 
requires an enormous number of RF components and bandpass filters, which leads to high cost. Besides, 
the number of the bands is fixed and the filter range is always preset. Thus the filter bank way is not 
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flexible. The other one is a wideband circuit using a single RF chain followed by high-speed digital 
signal processor (DSP) to flexibly search over multiple frequency bands concurrently |l9l. It is flexible to 
dynamic power spectrum density. High sampling rate requirement and the resulting large number of data 
for processing are the major problems IfTOl . 

Too high sampling rate requirement brings challenge to the analog-to-digital converter (ADC). And the 
resulting large amount of data requires large storage space and heavy computation burden of DSP. Since 
survey shows sparsity exists in the frequency domain for primary signal, compressive sensing (CS) can 
be used to effectively decrease the sampling rate [fTTI [fT^ fT^ . It assets that a signal can be recovered 
with a much fewer randomized samples than Nyquist sampling with high probability on condition that 
the signal has a sparse representation. 

In compressive wideband spectrum sensing (CWSS), analog-to-information converter (AIC) can be 
taken to obtain the random samples from analog signal in hardware as Fig. [T] shows lfT4l ifTSl . To get the 
spectrum estimation, there are mainly two groups of methods |fT3l . One group is convex relaxation, such 
as basis pursuit (BP) [fT6l IfTTll . Dantzig Selector (DS) [18J , and so on; the other is greedy algorithm, 
such as matching pursuit (MP) |fT91 , orthogonal matching pursuit (OMP) [EUll , and so on. Both of the 
convex programming and greedy algorithm have advantages and disadvantages when applied to difi'erent 
scenarios. A short assessment of their diff"erences would be that convex programming algorithm has a 
higher reconstruction accuracy while greedy algorithm has less computation complexity. In contrast to 
BP, basis pursuit denoising (BPDN) has better denoising performance [fTTI [flTTl . 

In this paper, the partial Fourier random samples are obtained via AIC with the measurement matrix 
generated by choosing part of separate rows randomly from the Fourier sampling matrix [fT4l . Based on 
the random samples, a generalized sparse constraint in the form of mixed '^2/'^! norm is proposed to 
enhance the recovery performance by exploiting the structure information. It encourages locally cluster 
distribution and globally sparse distribution. In the constraint, the estimated spectrum vector is divided 
into sections with different length according to the a priori information about fixed spectrum allocation. 



JOURNAL TITLE, VOL. X, NO. X, MONTH YEAR 4 

The sum of weighted ^2 norms of the sections is minimized. The weighting factor is iteratively updated 
as the reciprocal of the energy in the corresponding subband to get more democratical penalty of nonzero 
coefficients. Simulation results demonstrate that the proposed generalized sparse constraint based CWSS 
gets better performance than the traditional methods in spectrum reconstruction accuracy. 

In the rest of the paper, Section II gives the signal model; Section III states the classical CWSS methods. 
Section IV provides the generalized sparse constraint based CWSS methods; In section V, the performance 
enhancement of the proposed method is demonstrated by numerical experiments; Finally Section VI draws 
the conclusion. 

II. Signal Model 

According to the FCC report the allocated spectrum is in a very low utilization ratio. It means 
the spectrum is in sparse distribution. Recently a survey of a wide range of spectrum utilization across 6 
GHz of spectrum in some palaces of New York City demonstrated that the maximum utilization of the 
allocated spectrum is only 13.1%. It is also the reason that CR can work. Thus it is reasonable that only a 
small part of the constituent signals will be simultaneously active at a given location and a certain range 
of frequency band. The sparsity inherently exists in the wideband spectrum IfTOl Il22l ||23l Il24l ||25l ||26l 

112a m. 

An N X I signal vector x can be expanded in an orthogonal complete dictionary ^nxn, with the 
representation as 

Xa?x1 = ^WxwbiVxl (1) 

When most elements of the x 1 vector b are zeros, the signal x is sparse. When the number of nonzero 
elements of b is 5 (5 «c M < A''), the signal is said to be S'-sparse. 

In traditional Nyquist sampling, the time window for sensing is t e [0,Tq]. N samples are needed to 
recover the frequency spectrum r without aliasing, where Tq is the Nyquist sampling duration. A digital 
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receiver converts the continuous signal x(t) to a discrete complex sequence y, of length M. For illustration 
convenience, we formulate the sampling model in discrete setting as it does in IfTOl Il22l ||23l Il24l ||25l 
nasi [|21|I281: 

yt = AXf (2) 

where x, represents an A'^ x 1 vector with elements Xt[n] = x(t), t = hTq, n = I, - • • ,N, and A is an M x A'^ 
projection matrix. For example, when A = Fn with M = N, model ^ amounts to frequency domain 
sampling, where is the A'^-point unitary discrete Fourier transform (DFT) matrix. Given the sample 
set X, when M < N, compressive spectrum sensing can reconstruct the spectrum of r(t) with the reduced 
amount of sampling data. 

To monitor such a broad band, high sampling rate is needed. It is often very expensive. Besides, too 
many sampling measurements inevitably ask more storage devices and result in high computation burden 
for digital signal processors (DSP), while spectrum sensing should be fast and accurate. CS provides an 
alternative to the well-known Nyquist-Shannon sampling theory. It is a framework performing non-adaptive 
measurement of the informative part of the signal directly on condition that the signal is sparse ^T3\ . Since 
it is proved that x^ has a sparse representation in frequency domain. We can use an MxN random projection 
matrix to sample signals, i.e. = S^Xr, where M < N; Sc is a non-uniform subsampling or random 
subsampling matrix which is generated by choosing M separate rows randomly from the unit matrix I^y. 

The AIC can be used to sample the analog baseband signal x(t). One possible architecture can be 
based on a wideband pseudorandom demodulator and a low rate sampler [fT4l [fT5l . First we modulate 
the analogue signal by a pseudo-random maximal-length PN sequence. Then a low-pass filter follows. 
Finally, the signal is sampled at sub-Nyquist rate using a traditional ADC. It can be conceptually modeled 
as an ADC operating at Nyquist rate, followed by random discrete sampling operation [|T4ll . Then is 
obtained directly from continuous time signal x(t) by AIC. The details about AIC can be found in [[T4l 
[fT5l . Here we incorporate the AIC to the spectrum sensing architecture as Fig. [T] shows. 
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in. The Classical Compressive Wideband Spectrum Sensing 

CS theory asserts that, if a signal has a sparse representation in a certain space, one can use the 
random sampling to obtain the measurements and successfully reconstruct the signal with overwhelming 
probability by nonlinear algorithms, as stated in section II. The required random samples for recovery are 
far fewer than Nyquist sampling. 

To find the unoccupied spectrum for secondary access, the signal in the monitored band is down- 
converted to baseband. The analog baseband signal is sampled via the AIC that produces measurements 
at a rate below the Nyquist rate. 

Now we estimate the frequency response of x( t) from the measurement vector based on the transfor- 
mation equality y, = ScF^'r, where r is the x 1 frequency response vector (FRV) of signal x(t); is 
the NxN Fourier transform matrix; is the MxN matrix which is obtained by randomizing the column 
indices and getting the first M columns. 

Under the sparse spectrum assumption, the FRV can be recovered by solving the combinatorial opti- 
mization problem 

f = argmin ||r||o 

(3) 

s.t. (S,^F;^i)r = y, 

Since the optimization problem @ is nonconvex and generally impossible to solve, for its solution usually 
requires an intractable combinatorial search. As it does in [fTOll . BP is used to recover the signal: 

Tbp = argmin ||r||i 

(4) 

s.t. (sjF;^i)r = y, 

This problem is a second order cone program (SOCP) and can therefore be solved efficiently using standard 
software packages. 

BP finds the smallest norm of coefficients among all the decompositions that the signal is decomposed 
into a linear combination of dictionary elements (columns, atoms). It is a decomposition principle based 
on a true global optimization. 
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In practice noise exists in data. Another algorithm called BPDN has superior denoising performance 
than BP [21 J. It is a shrinkage and selection method for linear regression. It minimizes the sum of the 
absolute values of the coefficients, with a bound on the sum of squared errors. To get higher accuracy, 
we can formulate the BPDN based compressive wideband spectrum sensing (BPDN-CWSS) optimization 
model as: 



Tbpdn = argmm||r||i 

(5) 



s.t. ||(S,^F;^^)r-y, 

where 771 bounds the amount of noise in the data. The computation of the BPDN is a quadratic programming 



problem or more general convex optimization problem, and can be done by classical numerical analysis 
algorithms. The solution has been well investigated ETI ||29l [|30]| OTll . A number of convex optimization 
software, such as cvx [1321 , SeDuMi [i33il and Yalmip [13?ll , can be used to solve the problem. 

IV. The Proposed Compressive Wideband Spectrum Sensing 

Among the classical sparse signal recovery algorithms, BPDN achieves the highest recovery accuracy 
[fT3l . However, it only takes advantage of sparsity. In wideband CR application, additional a priori 
information about the spectrum structure can be obtained. The further exploitation of structure information 
would give birth to recovery accuracy enhancement [|28l [|35l [|36l . Besides, It is well-known that the 
minimization of norm is the best candidate for sparse constraint. But in order to reach a convex 
programming, the "^o norm is relaxed to ^1 norm, which leads to the performance degeneration [|37l . 
Here a weighting formulation is designed to democratically penalize the elements. It suggests that large 
weights could be used to discourage nonzero entries in the recovered FRY, while small weights could be 
used to encourage nonzero entries. To get the weighted values, a simple iterative algorithm is proposed. 

A. Wideband spectrum sensing for fixed spectrum allocation 

The classical algorithms reconstruct the commonly sparse signal. However, in the coarse wideband 
spectrum sensing, the boundaries between different kinds of primary users are fixed due to the static 
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frequency allocation of primary radios. For example, the bands 1710 - 1755 MHz and 1805 - 1850 MHz 
are allocated to GSM 1800. Previous CWSS algorithms did not take advantage of the information of fixed 
frequency allocation boundaries. Besides, according to the practical measurement, though the spectrum 
vector is sparse globally, in some certain allocated frequency sections, they are not always sparse. For 
example, in a certain time and area, the frequency sections 1626.5 - 1646.5 MHz and 1525.0 - 1545.0 
MHz allocated to international maritime satellite are not used, but the frequency sections allocated to 
GSM 1800 are fully occupied. The wideband FRY is not only sparse, but also in sparse cluster distribution 
with different length of clusters. It is the generalization of the so called block-sparsity [|35l ||36l . This 
feature is extremely vivid in the situation that most of the monitored primary signals are spread spectrum 
signals. 

Previous classical CWSS does not assume any additional structure on the unknown sparse signal. 
However in the practical application, the signal may have other structures. Incorporating additional 
structure information would improve the recoverability potentially. 

Block-sparse signal is the one whose nonzero entries are contained within several clusters. To exploit 
the block structure of ideally block-sparse signals, ^2/^1 optimization was proposed. The standard block 
sparse constraint (SBSC) in the form of '^2/'^! optimization can be formulated as ll35l [f36ll : 



where K is the number of the divided subbands; do is the length of the divided blocks. Extensive 
performance evaluations and simulations have demonstrated that as do grows the algorithm significantly 
outperforms standard BP algorithm [|36ll . 

However, in the standard ^2/^1 optimization, the estimated sparse signal is divided with the same block 
length, which mismatches the practical situation that the values of the length of the spectrum subbands 
allocated to different radios can not be all the same. Besides, the constraint in Q does not incorporate 
the denoising function. 




) 



(6) 



JOURNAL TITLE, VOL. X, NO. X, MONTH YEAR 9 

To further enhance the performance of CWSS, the fixed spectrum allocation information can be in- 
corporated in the CWSS algorithm. Based on the a priori information about boundaries, the estimating 
PSD vector is divided into sections with their edges in accordance with the boundaries of different types 
of primary users by fixed spectrum allocation. In the BPDN-CWSS, the minimization of the standard 
"^i-norm constraint on the whole FRV is replaced by the minimization of the sum of the ^2 norm of 
each divided section of the FRV to encourage the sparse distribution globally while blocked distribution 
locally. As it combines ^\ norm and ^2 norm to enforce the sparse blocks with different block lengths, the 
new CWSS model, in the name of variable-length-block- sparse constraint based compressive wideband 
spectrum sensing (VLBS-CWSS), can be formulated as: 



min(||ri||2 + ||r2||2 + ---||r^|l2) 

(7) 

s. t. ||y, - S,F^ir||2 < 772 

where ri, r2, ... , r^^ are K sub- vectors of r corresponding to di, di, ... , dK~i which are the boundaries 
of the divided sections. 772 bounds the amount of noise in the data. It can be formulated as: 



fl ■■■ fell ''' ^dK-\ + l 



(8) 



V ri Tk 

Since the objective function in the VLBS-CWSS © is convex and the other constraint is an affine, it 
is a convex optimization problem. It can also be solved by a host of numerical methods in polynomial 
time. Similar to the solution of the BPDN-CWSS ©, the optimal r of the VLBS-CWSS dV]) can also be 
obtained efficiently using some convex programming software packages. Such as cvx [|32]| . SeDuMi [|33]| . 
and Yalmip BMl, etc. 

After we get r from ([8]), power spectrum can be obtained. Several ways can indicate the spectrum 
holes, such as energy detection [[27ll . edge detection [fTOll . and so on. For example, in energy detection we 
will calculate \\rk\\2, k = \, 2, ... , K. Comparing it with an experimental threshold, the spectrum holes for 
dynamic access can be clearly given. The energy detection will be used in numerical simulations. 
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B. Enhanced variable-length-block-sparse spectrum sensing 

In sparse constraint, %) norm minimization is relaxed to norm at the cost of bringing the dependence 
on the magnitude of the estimated vector. In the norm minimization, larger entries are penalized more 
heavily than smaller ones, unlike the more democratic penalization of the norm. Here in the the VLBS 
constraint, to encourage sparse distribution of the spectrum in the global perspective, the norm of a 
series of the norm is minimized. Similarly, the dependence on the power in each subband exits. 

To deal with this imbalance, the minimization of the weighted sum of the ^2 norm of each blocks 
is designed to more democratically penalize. The new weighted VLBS constraint based compressive 
wideband spectrum sensing (WVLBS-CWSS) can be formulated as: 



min (wi ||ri||2 + W2 \\r2\\2 + ■ ■ ■ + w/f llr^f ly 

(9) 



s.t. y, - S,F.ir I < 773 



N * II2 

T 

■ Wi 



Wi W2 ■ ■ ■ Wk 



where ri, r2, ... , Yk are defined as (|8]); 773 bounds the amount of noise; w 
depends on > O,for / = I, - - ,K, where corresponds to the power of the primary user exists in the 
i-th subband. 

Obviously, the object function of the WVLBS-CWSS Q is convex. It is a convex optimization problem. 
In principle this problem is solvable in polynomial time. 

To realize the WVLBS-CWSS (|9l), the weighting vector w should be provided. As it is defined before, 
the computation of the weight w, is in fact the computation of the p,. Here a practical way to iteratively 
set the Pi is proposed. At each iteration, the pi is the sum of the absolute value of frequency spectrum 
vector in the corresponding subband. It can be formulated as: 

Pui = h-u\l 

111 ^^Q^ 
= \ft-l, £/,_i + l| + • • • + di\ 

where r,_i , is the i-th sub-vector as in ([8]) at the {t-\)-th iteration; r,_i • • • , ^. are the elements 

of the sub-vector r^.i After getting the pi, the weighting vector w can be formulated. Here we can get 
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it by 

Wi = -^ (11) 

Pi + 6 

where a small parameter 5 > in (fTTI) is introduced to provide stability and to ensure that a zero-valued 
component in pi does not strictly prohibit a nonzero estimate at the next step. 

The initial condition of the recursive relation is w, = 1, for all i = 1, ...,K. That means in the first step, 
all the blocks are weighted equally. Along with the increase of the iteration times, larger values of p, are 
penalized lighter in the WVLBS-CWSS ^ than smaller values of pi. To terminate the iteration at the 
proper time, the stopping rule can be formulated as 



||r, -r,_i||2 <£ (12) 

where is the estimated FRV at the t-th iteration; s bounds the iteration residual. 

The initial state of the iterative algorithm is the same with the VLBS-CWSS Q. To make a difference. 
The iterative reweighted algorithm is named as enhanced variable-length-block-sparse constraint based 
compressive wideband spectrum sensing (EVLBS-CWSS). 

V. Simulation Results 

Numerical experiments are presented to illustrate performance improvement of the proposed EVLBS- 
CWSS for CR. Here we consider a base band signal with its frequency range from Hz to 500 MHz as 
Fig. [2] shows. The primary signals with random phase are contaminated by a zero-mean additive white 
Gaussian noise (AWGN) which makes the signal to noise ratio (SNR) be 11.5 dB. Four primary signals 
are located at 30 MHz - 60 MHz, 120 MHz - 170 MHz, 300 MHz - 350 MHz, 420 MHz - 450 MHz. 
Their corresponding frequency spectrum levels fluctuate in the range of 0.0023 - 0.0066, 0.0016 - 0.0063, 
0.0017 - 0.0063, and 0.0032 - 0.0064, as Fig. [3] shows. Here we take the noisy signal as the received 
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signal x(t). As CS theory suggests, we sample x(t) randomly at the subsampling ratio 0.40 via AIC as 
Fig. [TJ The resulted sub-sample vector is denoted as 

To make contrast, with the same number of samples. The amplitude of frequency spectrum estimated 
by different methods are given in Fig. |4l FigjS] and Fig. [6l Fig. |4] shows the result estimated by the 
standard BPDN-CWSS ^ where rji is chosen to be 0.1 \\yt\\2 with 1000 tries averaged; Fig. [5] does it by 
the VLBS-CWSS © where j]2 is chosen to be 0.2 Hy^Hj; Fig. [6] does it by the proposed EVLBS-CWSS 
(|9]) where 773 is chosen to be 0.2 WytWj, and 6 is chosen to be 0.001. 

Fig. |6] shows that the proposed EVLBS-CWSS gives the best reconstruction performance. It shows 
that there are too many fake spectrum points in the subbands with no active primary signal in Fig |4] 
which is given by the standard BPDN. The noise levels of the spectrum estimated by the B-CWSS and 
the VLBS-CWSS are high along the whole monitored band. For the VLBS-CWSS, as in Fig. [51 it has 
considerable performance improvement, but the noise level in part of the inactive subbands is still high. 
Some of the estimated spectrum in the inactive subband is a little too high. However, in Fig. [6l the four 
occupied bands clearly show up; the noise levels in the inactive bands are quite low; the variation of the 
spectrum levels in the boundaries of estimated spectrum are quite abrupt and correctly in accordance with 
the generated sparse spectrum in Fig. [2l which would enhance the edge detection performance much. 
Therefore, the proposed EVLBS-CWSS outperforms the standard BPDN-CWSS and the VLBS-CWSS 
for wideband spectrum sensing. 

Apart from the edge detection, energy detection is the most popular spectrum sensing approach for CR. 
To test the CWSS performance by energy detection, 1000 Monte Carlo simulations are done with the same 
parameters above to give the results of average energy in each section of the divided spectrum vector with 
the BPDN-CWSS ©, the VLBS-CWSS O and the EVLBS-CWSS @. The parameter setting is same 
as before. The simulated monitored band is divided into 9 sections as Fig [2l The total energy with each 
CWSS method is normalized. Table |I] presents the average energy in each subband with different recovery 
methods, when there are 4 active bands and the sub-sampling ratio is 0.40; Table HI] does when there are 
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3 active bands and the sub-sampling ratio is 0.40; Table [nl] does when there are 3 active bands and the 
sub-sampling ratio is 0.35; Table |IV] when there are 2 active bands and the sub-sampling ratio is 0.30. 
For the EVLBS-CWSS, it is obvious that the estimated noise energy of inactive bands is much smaller 
that the other two. To quantify the performance gain of EVLBS-CWSS against others, after normalizing 
the total energy of the spectrum vectors, we define the energy detection performance enhancement ratios 
(EDPER) of VLBS-CWSS and EVLBS-CWSS against BPDN-CWSS for the k-th subband as: 



R\{k) 



|j.VLBS |H_||j.S/>DA'||- 

* — for active subbands 



W-BPDN W- UfVLBS 
I't lb W k Il2 

IF* lb 



(13) 

, for inactive subbands 



R2 (k) = 



\j.EVLBS U'^Uj.BPDNU- 

— ^ * — —, for active subbands 



I Ft lb 

J.BPDN \ U^EVLBS U'^ 

- — ^ —, for inactive subbands 



(14) 



\-BPDN\ 

I'a I 



where r^^^^^, r^^^-^ and r^™^ represent values of estimated frequency spectrum vectors in the k-th 
subband via EVLBS-CWSS, VLBS-CWSS and BPDN-CWSS, respectively. These performance functions 
can quantify how much energy increased to enhance the probability of correct energy detection of the 
active primary bands and how much denoising performance is enhanced. The values of EDPER in Table 
i Table HIl Table |lll] and Table |IV] clearly tell the improvement of the proposed EVLBS-CWSS against 
VLBS-CWSS and BPDN-CWSS methods. 

To further evaluate the performance of EVLBS-CWSS, when the number of active bands is 4 and sub- 
sampling ratio is 0.40, the residuals ||r,-rf_i||2 for 1000 Monte Carlo simulations are measured. Using 
the unnormalized received signal, the measured average power of the random samples is 29533. From 
? = 2 to ? = 8, the residuals are 361.5066, 261.6972, 55.0035, 17.9325, 15.0799, 13.4075 and 12.6189. It 
shows the iteration is almost convergent at t = 5. The iteration would bring the increase of computation 
complexity, but the performance enhancement is obvious and worthwhile. 

The enhancement of spectrum estimation accuracy qualifies the proposed EVLBS-CWSS as an excellent 
candidate for CWSS. 
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VI. Conclusion 

In this paper, CS is used to deal with the too high sampling rate requirement problem in the wideband 
spectrum sensing for CR. The sub-Nyquist random samples is obtained via the AIC with the partial 
Fourier random measurement matrix. Based on the random samples, incorporating the a priori information 
of the fixed spectrum allocation, an improved sparse constraint with diff"erent block length is used to 
enforce locally block distribution and globally sparse distribution of the estimated spectrum. The new 
constraint matches the practical spectrum better. Furthermore, the iterative reweighting is used to alleviate 
the performance degeneration when the norm minimization is relaxed to the ^2/% one. Because 

the a priori information about boundaries of diff'erent types of primary users is added and iteration is 
used to enhance the VLBS constraint performance, the proposed EVLBS-CWSS outperforms previous 
CWSS methods. Numerical simulations demonstrate that the EVLBS-CWSS has higher spectrum sensing 
accuracy, better denoising performance, etc. 
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Fig. 1. The proposed compressive wideband spectrum sensing structure. 
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Fig. 2. Thie normalized spectrum of noiseless active primary signals in the monitoring band. 
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Fig. 3. The normalized spectrum of noisy active primary signals in the monitoring band. 
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Fig. 4. The compressive wideband spectrum estimation via BPDN-CWSS. 
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Fig. 5. The compressive wideband spectrum estimation via VLBS-CWSS. 
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TABLE I 

The total energy in each subband with the three CWSS methods and the values of EDPER, when there are 4 active bands and the 

sub-sampling ratio is 0.40. 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Real PSD 





0.4747 





5303 





0.5107 





0.4823 





BPDN-CWSS 


0.1149 


0.3820 


0.1752 


0.4734 


0.3184 


0.4780 


0.2333 


0.4026 


0.1994 


VLBS-CWSS 


0.0000 


0.2447 


0.0000 


0.5101 


0.4220 


0.5833 


0.0005 


0.4020 


0.0000 


EVLBS-CWSS 


0.0000 


0.2681 


0.0000 


0.5396 


0.1897 


0.6361 


0.0000 


0.4431 


0.0000 


Rl 


100 % 


-35.94 % 


100 % 


7.75 % 


-32.54 % 


22.03 % 


100 % 


0.15 % 


100 % 


R2 


100 % 


-29.82 % 


100 % 


13.98 % 


40.42 % 


33.08 % 


100 % 


10.06 % 


100 % 
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TABLE n 

The total energy in each subband with the three CWSS methods and the values of EDPER, when there are 3 active bands and the 

sub-sampling ratio is 0.40. 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Real PSD 


0000 


0000 


0000 


0.5998 


0000 


0.6171 


0000 


0.5093 


0000 


BPDN-CWSS 


0.2489 


0.1221 


0.1704 


0.5080 


0.2526 


0.5676 


0.1637 


0.4867 


0.1741 


VLBS-CWSS 


0.0000 


0.0000 


0.0000 


0.5642 


0.2544 


0.6806 


0.0000 


0.3922 


0.0000 


EVLBS-CWSS 


0.0000 


0.0000 


0.0000 


0.6029 


0.0027 


0.5951 


0.0000 


0.5313 


0.0000 


Rl 


100 % 


100 % 


100 % 


11.06 % 


-0.71 % 


19.91 % 


100 % 


-19.42 % 


100 % 


R2 


100 % 


100 % 


100 % 


18.68 % 


98.93 % 


4.84 % 


100 % 


9.16 % 


100 % 
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TABLE m 

The total energy in each subband with the three CWSS methods and the values of EDPER, when there are 3 active bands and the 

sub-sampling ratio is 0.35. 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Real PSD 


0000 


0000 

\J t\J \J\J \J 


0000 


0.5997 


0000 


0.6171 


0000 


0.5094 


0000 


BPDN-CWSS 


0.1966 


0.1766 


0.2002 


0.5912 


0.3480 


0.4420 


0.2573 


0.3482 


0.1915 


VLBS-CWSS 


0.0000 


0.0000 


0.0000 


0.8035 


0.4190 


0.3593 


0.0000 


0.2231 


0.0000 


EVLBS-CWSS 


0.0000 


0.0000 


0.0000 


0.7572 


0.1258 


0.5095 


0.0000 


0.3887 


0.0000 


Rl 


100 % 


100 % 


100 % 


35.91 % 


-20.40 % 


-18.71 % 


100 % 


-35.93 % 


100 % 


R2 


100 % 


100 % 


100 % 


28.08 % 


63.85 % 


15.27 % 


100 % 


11.63 % 


100 % 
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TABLE rV 

The total energy in each subband with the three CWSS methods and the values of EDPER, when there are 2 active bands and the 

sub-sampling ratio is 0.30. 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Real PSD 


0000 


0000 


0000 


0.7918 


0000 


0000 


0000 


0.6107 


0000 


BPDN-CWSS 


0.1710 


0.0736 


0.1733 


0.5836 




0.1907 


0.2341 


0.6021 


0.2565 


VLBS-CWSS 


0.0000 


0.0000 


0.0000 


0.7768 


0.2151 


0.0000 


0.0000 


0.5907 


0.0000 


EVLBS-CWSS 


0.0000 


0.0000 


0.0000 


0.7697 


0.0012 


0.0000 


0.0000 


0.6387 


0.0000 


Rl 


100 % 


100 % 


100 % 


33.10 % 


21.55 % 


100 % 


100 % 


-1.89 % 


100 % 


R2 


100 % 


100 % 


100 % 


19.60 % 


74.93 % 


100 % 


100 % 


6.08 % 


100 % 



